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ABSTRACT 

We compute photoionization models for the giant extragalactic H ll region 
NGC 5461, and compare their predictions to several observational constraints. 
Since we aim at reproducing not only the global properties of the region, but 
its local structure also, the models are constrained to reproduce the observed 
density profile, and our analysis takes into consideration the bias introduced by 
the shapes and sizes of the slits used by different observers. We find that an 
asymmetric nebula with a gaussian density distribution, powered by a young 
burst of 3.1 Myr, satisfactorily reproduces most of the constraints, and that the 
star- formation efficiency inferred from the model agrees with current estimates. 
Our results strongly depend on the assumed density law, since constant density 
models overestimate the hardness of the ionizing field, affecting the deduced 
properties of the central stellar cluster. We illustrate the features of our best 
model, and discuss the possible sources of errors and uncertainties affecting the 
outcome of this type of studies. 

Subject headings: H ll regions — ISM: abundances — ISM: individual (MlOl, 
NGC5461) — Stars: formation. 



1. INTRODUCTION 

Tailored photoionization models have proven to be a useful tool for the understanding of 
star- formation (SF) regions (e.g., Garcia- Vargas et al. 1997; Gonzalez-Delgado & Perez 2000; 
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Luridiana, Peimbert, & Leitherer 1999; Luridiana 1999; Stasinska & Schaerer 1999). This 
has been especially true since spectral energy distributions (SEDs) for population synthesis 
models have become available, providing a much better approximation to real ionizing spectra 
than the naive one-spectrum models. 

Many uncertainties still affect the output of photoionization models: incomplete stellar 
tracks, unknown geometry both of the stellar source and the ionized nebula, processes -other 
than photoionization- participating in the gas thermal balance, to mention only a few of 
them. An excellent review on this kind of uncertainties can be found in Stasinska (2000). 
An additional source of uncertainty lies in the stochasticity of the Initial Mass Function 
(IMF) in real clusters, which leads to fluctuations of the number of stars of a given mass 
value around the average (analytical) value (Cervino, Luridiana & Castander 2000). This 
source of uncertainty is intrinsic to population synthesis models and cannot be removed, 
and it especially affects low-mass star clusters as it depends on statistics; it is probably not 
relevant for the case of NGC 5461, as it will be shown in Section 5.2. 

In spite of these uncertainties, the computation of tailored photoionization models can 
provide many insights into the physical processes going on in photoionized regions. In the 
present work we describe a selected photoionization model of NGC 5461, and illustrate the 
properties that can be inferred for the region, emphasizing the uncertainties involved and the 
problems still unsolved. The predictions of the model are compared to a large set of obser- 
vational constraints, including both global (e.g., the total emitted H/? flux) and more local 
properties (e.g., the line intensity ratios for diflcrcnt apertures). The comparison between 
the predictions and the observations is always made after correcting the model output for the 
size and shape of the aperture used in each observation. Furthermore, we did not enforce any 
a priori density law, but rather determined a gas distribution yielding self-consistent pre- 
dictions in agreement with resolved radial observations of the /(A6717)//(A6731) doublet. 
This is by far the most innovative point in our modeling procedure, and, quite certainly, one 
of the most important: in fact, in the following we will demonstrate that our results strongly 
depend on the assumed density law, and that simpler, not-tailored gas distributions, such as 
constant density models taken from large photoionization grids, overestimate the hardness 
of the radiation field, leading to a significant bias in the deduction of the properties of the 
ionization source. 



2. GENERAL PROPERTIES 



NGC 5461 is a giant cxtragalactic H ll region (GEHR), located in one arm of the spiral 
galaxy M 101 (NGC 5457). NGC 5461 has been the object of several studies (Israel, Goss, & 
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AUen 1975; Sandage & Tammann 1976; Rayo, Peimbert, & Torres-Peimbert 1982; McCall, 
Rybski, & Shields 1985; Evans 1986; Melnick ct al. 1987; Skillmann & Israel 1988; Torres- 
Pcimbert, Peimbert, & Fierro 1989; Castaiieda, Vilchcz, & Copetti 1992; Rosa & Bcnvenuti 
1994; Williams & Chu 1995; Kennicutt & Garnctt 1996; Garnett et al. 1999). In the 
modeling, we used as many data as possible from these sources, unless relevant information 
from the reference paper was missing. Additionally, we used high-resolution spectroscopic 
data taken in June 1996 with the 2.1 m telescope at the Observatorio Astronomico Nacional 
de S. Pedro Martir (Luridiana, Esteban, & Peimbert, in preparation). In the following 
sections we will discuss the determinations of relevant parameters made by different authors. 

2.1. Radius 

The radius of an extended object is a somewhat tricky parameter, since it depends on 
the lower limit set on the observed flux, and on the frequency considered. Castaiieda et 
al. (1992) reproduce the Ho; brightness profile of NGC 5461 with the superposition of two 
gaussian density distributions, of characteristic radius ro = 1.3", separated by 2" in the plane 
of the sky. Williams & Chu (1995) show an uncalibrated Ha image, in which the maximum 
emission comes from an elongated region of about 30" x 15". Kennicutt & Garnett (1996), 
based on data by Scowen, report a diameter of 30" in Ha. Israel et al. (1975) estimate a 
size of 66" x 25" in Ha, and find that the radio emission follows the same distribution; they 
state that more than 50% of the flux is emitted in a bulge with a FWHM of 5". 

We took into account these data to have a gross idea of the radius of the region, which 
we estimate to be about 14". At a distance of 7.4 Mpc (Sandage & Tammann 1976, with the 
suggested correction by dc Vaucouleurs 1978) this corresponds to a linear radius of 500 pc. 
In all our calculations, we always ensured that the total radius was not too different from 
this reference value. 



2.2. Electron Density and Filling Factor 

Castaiieda et al. (1992) found that a model following the density law A^c = A'^e'^*^''/''")^ + 
A'te^*-*-''^^-'/''"-*^ reproduces the observed Ha brightness and /(A6717)//(A6731) profiles. Our 
models are constrained to reproduce their observed /(A6717)//(A6731) profile, however we 
did not stick to the density law they proposed for the following reasons: 

1. The density law cannot be calculated a priori knowing only the observed /(A 6717) / /(A 6731) 
profile, but it must rather be determined self-consistently with the overall ionization 
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structure. Elliott & Meaburn (1973) give an approximate expression to estimate the 
[O ll] density of a dust-free nebula: 

iV.(OIIl) = f|=<^. (1) 

where V is the volume and £"3729 ("i^), -E'3726('^) represent the local emissivities of [O II] 
A 3729 and [O ll] A 3726. Castaneda et al. (1992) generalize this expression to the 
[S 11] case; however, they neglect to include the S ll abundance in the integrations, 
which does not cancel out since it strongly varies from point to point; in this respect, 
it is sufficient to note that the presence of O ill in the spectrum implies the presence 
of S IV, which has roughly the same ionization potential. The fact that the S ll 
abundance is not spatially constant is also noted by Castaneda et al. (1992), when 
they observe that S ll is not the most important ionization state of sulfur in a GEHR, 
and that it is present only in the outskirts of the nebula. 

Stated otherwise, the observed /(A6717)//(A6731) profile at a given point does not 
only depend on the physical conditions (electron density and temperature), but also on 
the ionization conditions found at each projected point: e.g., the [S II] emission from 
a central high-density zone can be diluted off by the contribution of an extended low- 
density halo, being invisible in projection; or, if the highest-density zone is highly ion- 
ized, it does not contain any S ll and gives no contribution to the /(A 6717)//(A 6731) 
ratio. 

2. For simplicity, we chose a centrally symmetric density law. This choice implies cither 
neglecting the secondary brightness knot, seen a.t d = —2" in Fig. 13 of Castaneda et 
al. (1992), or adopting an approximate scheme to model the asymmetry. We adopted 
the second solution, as it will be explained in the following. 



Following Castaneda et al. (1992), we adopted a gaussian form for the density law. 
Within this choice, two alternatives are possible: a) Either the ionizing source is located at 
the brightest spots, or b) it is shifted with respect to them. In principle, both solutions are 
possible, since the emitted intensity depends both on the density and on the number of ion- 
izing photons striking a given point. However, the /(A6717)//(A6731) ratio profile indicates 
that the source is most probably located somewhere between the two emission maxima. In 
fact, in the innermost, high-ionization zone of the H II region no S II is expected (hence no 
[S 11] emission), so that a local minimum in density (local maximum in /(A 671 7) //(A 6731)) 
should be observed. Instead, the [S ll] ratio shows a definite absolute minimum at the 
brightest knot location, and a secondary local minumum at the secondary knot. This prob- 
ably implies that the stellar cluster is located between the two, and has blown away and/or 
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consumed out the neighboring gas, leaving a density depression surrounded by the density 
peaks seen in projection. A continuum image of the region, showing the position of the stars 
with respect to the gas, would be of great help in confirming this point. 

Based on these considerations, we modeled the region with off-center gaussians, following 
the density law 



where /(r) = iVoe-«'-^)/'"o)', Nf"" = 50 cm'^, and 1.0" < S < 2.0" and No are tuned to 
meet the observational constraints. Since the brightness knots produced in projection by a 
model of this kind are identical, whereas the peak intensity of the secondary knot in NGC 
5461 is about one half that of the primary one (Fig. 13 in Castaiieda et al. 1992), we 
calculated asymmetric models, by adding the two halves of two different models, ionized by 
the same source but differing in the density distributions. This configuration neglects the 
pressure gradients at the border between the two halves and assumes that the diffuse fields 
are directed outwards. 

The filling factor can be estimated by means of the relation: 



where N^{rms) is the root mean square electron density, and N^{FL) is the electron density 
determined through a forbidden-line ratio. With the values given by Torres-Peimbert et al. 
(1989) for A^'g^(rms) and A^|([S ll]), one finds e = 0.004; using more recent atomic data for 
[S ll] (Pradhan & Peng 1995; Keenan et al. 1996) we find e = 0.002 (see Table 3). 

In our numerical experiments, we found it difficult to simultaneously obtain the desired 
brightness and 7(A 6717)//(A 6731) profiles with two different A^o values for the two halves: 
instead, we found rather good fits by adopting equal A^o values and different filling factors 
(e) on the two sides. Fig. 1 sketches a section of a model of this kind, with the dot size pro- 
portional to the density, and the dot spacing proportional to the fiUing factor. The assumed 
positions of some of the slits used to observe the region are also shown (see Sections 2.5 
and 3.2). 

As for the filling factor behavior in each half, we ignore whether the choice of keeping it 
constant is physically based, nor whether it is the only possible solution: rather, we chose it 
as the simplest hypothesis leading to self-consistent predictions. A thorough discussion on 
the physical meaning of such configuration is out of the scopes of the present paper. 




for /(r) > Nf 
for /(r) < N^' 



(2) 



N^{rms)^eN^{FL), 



(3) 
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2.3. g(HO) 

The rate of ionizing photons emitted by the ionizing stars can be computed by means 
of the expression: 

^ ^ hTv 

where 7*°*(H/9) stands for the total dereddened H/3 intensity, measured in erg sec~^cm~^, d 
is the distance to the region, q;b is the total recombination coefficient to all levels but the 
first for hydrogen case B, and q;(H/3) is the effective recombination coefficient for H/3. The 

assumptions underlying this expression arc that no photons leak out nor arc absorbed by 
dust, and that case B holds. The effects of relaxing the first assumption depend on the 
constraints adopted on the models (see Section 6.3). 

(5(H°) values determined through observations are generally underestimated, since most 
slits sample only a part of the nebula, so that the observed values of /(H/5) are smaller 
than P"^{Hf3). To overcome this uncertainty, we considered as many data as possible, in a 
sequence of larger and larger slits. The first column of Table 1 lists the logarithmic values 
of the observed intensity of H/3 before reddening correction reported by different authors, 
in units erg sec~^cm~^. The second column shows the logarithmic reddening correction 
C(H/3), and the third column lists the corresponding (5(H°), homogeneously calculated by 
us adopting d = 7.4 Mpc and = 9, 200 K. 

We adopted for our models the maximum observed value Log Q(H°) = 52.54. Inciden- 
tally, this table illustrates the risks of modeling a region using only one set of observational 
constraints: since the local ionization degree directly depends on the rate of ionizing photons 
emitted, models in which the number of ionizing photons emitted is substantially underesti- 
mated cannot be realistic. 



2.4. Observed Equivalent Width of H/3 

The fourth column of Table 1 lists some of the values of the equivalent width of H/3 
{EW{}iP)) found in the literature. Though EW{}iP) is a relevant quantity, being at the 
same time significative and easy to measure, some caution should be taken before relating 
it to the age of the source. First of all, there could be a continuum emission contribution 
from older stars, not related to the present burst of star formation, lowering EW{}i(5) and 
thus mimicking an older age for the burst. Second, the equivalent width is, ideally, a global 
parameter, while the observed value is local since the slit samples only a part of the nebula. 
If the stars are concentrated in the center of the region, we should expect that the bigger 



Table 1. Compiled Parameters of NGC 5461 and Slits' Features. 



Measured Quantities Slits' Features 



Reference Log F(H/?) C(H/3) Log Q(HO)'* EW{Il/3) Shape Size Label 





(erg s ^) 


(dex) 


(photons s ^) 


(A-i) 




(") 


(pc)- 




Rayo et al. 1982 


-12.52 


0.75 




195 


Rectangular 


3.8 X 12.4 


137 X 446 




McCall et al 1985 




(0.60)*^ 




249 


Rectangular 


Several slits 




Mclnick ct al. 1987 


(-11.99)'= 


0.30 


52.44 


223 


Circular 


Unknown*^ 




Skillmann & Israel 1988 


-12.16 


0.57 


52.54 




Circular 


20 


© 720 


E 


Torres-Peimbert et al. 1989 


-12.46 


0.60 


52.27 


195 


Rectangular 


3.8 X 12.4 


137 X 446 


B 


Kennicutt & Garnctt 1996 




0.36 




239 


Circiilar 


(0 10.3?)« 


(0 371?) 




Luridiana et al. in preparation^ 


-12.39 


0.43 


52.18 


157 


Rectangular 


4.0 X 13.3 


144 X 479 


C 


Luridiana et al. in preparation^ 










Rectangular 


4.0 X 39.9 


144 X 1436 


D 


Garnett et al. 1999 




(0.57)'' 






Rectangular 


0.86 X 0.86 


31 X 31 


A 



''Calculated for an assumed distance of 7.4 Mpc. 

^Derived from the published F(Ha)/F(H/3), adopting the extinction law by Whitford (1958). 
^Derived from the published .F(H/3) = ATr(fF{E(3), and d = 6.9 Mpc. 

'^Greater than the halo diameter, but not specified by the authors. 
®Not clearly specified in the paper. 
^Bluc side of the spectrum (3300-6600 A). 
sRed side of the spectrum (6300-9100 A). 
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the slit, the higher the observed EW{}iP). This trend is just the opposite of what is found 
in Table 1, a not unusual finding in the study of GEHRs (see, e.g. the case of NGC 2363: 
Gonzalez-Delgado et al. 1994), probably implying continuum contamination by sparse low- 
mass stars. A third problem is related to the geometry of the source. For a given stellar 
population, the observed EW{B.f3) value is smaller for density bounded (as compared to 
radiation bounded) nebulae, since the /(H/3) value superposed to a given continuum is lower 
in the former case (this is probably one of the reasons underlying the apparent lack of very 
young starburst with EW(liP) values around 1000). 

With these uncertainties in mind, wc can assume the maximum observed value, £'M^(H/?)=249 
A~^, as a representative value for the stellar population responsible of the ionization, corre- 
sponding to ages of order 2.5 Myr ^ t ^ 3.5 Myr for the instantaneous burst SF (IBSF) law. 
This age range agrees with the values determined by Rosa & Benvenuti (1994) adopting a 
Miller-Scalo IMF. In the continuous SF (CSF) case, the mentioned EW{}if3) value would 
correspond to the age range 4.5 Myr ^ t ^ 5.5 Myr. 

2.5. Apertures 

The sixth to eighth columns of Table 1 hst the features of the sht apertures used by 
each observer, with the linear dimensions homogeneously calculated for an assumed distance 
of 7.4 Mpc to the parent galaxy. The last column hsts the labels used throughout the paper 
to refer to some of these slits. 



2.6. Chemical Abundances 

The adopted chemical abundances are hsted in Table 2. The abundance values for He, 
O, N, Ne, S and Ar used are those derived by Torres- Peimbert et al. (1989) through the 
empirical electron-temperature based method, since the previously published data (Rayo et 
al. 1982; McCall et al. 1985) were not corrected for the non-hnearity of the detectors (see 
Peimbert & Torres- Peimbert 1987). 

The C abundance is calculated averaging the two C/0 ratios (-0.03 and -0.37) by Garnett 
et al. (1999), corresponding to the cases Ry = Ay/E{B — V)— 3.1 and 5.0, and assuming 
log 0/H = -3.61 as in Torres-Peimbert et al. (1989). 

The total heavy element abundance is Zg^s — 0.0066 for the t'^ — case, where t'^ is 
the temperature- fluctuation parameter (Peimbert 1967). In our numerical experiments we 
also explored higher metallicity values, but found it difficult to discriminate among them 
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for reasons which will be explained in the following. Thus, the standard t'^ — value was 
adopted in the selected model presented in this paper. 



The photoionization model presented in Section 5 has been calculated assuming Zgas = 
0.0066 and = 0.008. The = 0.008 value assumes that approximately 20% of the 
heavy elements are locked in dust in the gas (Esteban et al. 1998); it also has the (purely 
technical) virtue of allowing us to avoid making interpolations between spectra of different 
metallicities, since it is one of the five Z values of the Starburst99 library. Nevertheless, a 
higher metallicity is in principle possible, given the following considerations: 

a) If there are temperature fluctuations in the object, the metallicity is underestimated 
if these are not taken into account (Peimbert 1967). 

b) The metallicity can be estimated from the Log R23 vs. [0/H] diagram (Pagel et al. 
1979; McGaugh 1991) when the ionization parameter U is known. With the oxygen line 
intensities listed in Table 3 one obtains 0.78 < Logi?23 < 0.84. Using the relationship 
by Diaz et al. (1991) between U and the intensity of the sulfur lines: 



one obtains 0.004 <U< 0.019, with the exact value depending on the atomic param- 
eters (Pradhan & Peng 1995; Shaw & Dufour 1994) and the data set used. With these 
R23 and U values, the calibration by McGaugh (1991) gives -3.4 < [0/H] < -3.2, 
corresponding to the range 0.010 < Zg^s < 0.017. 

c) Our model adjusts the R23, implying Zgg^ ~ 0.0066 in apparent contradiction with 
McGaugh's cahbration. However, it should be taken into account that NGC 5461 hes 
close to the bend of the it!23 diagram, so that the Z solutions are nearly degenerate 
with respect to i?23- Furthermore, C/ is a parameter which should be used with caution: 
though it is usually treated as a global property of the region, it is rather a local 
quantity. The observed U values are situated somewhere in between a local and a 
global definition, since they are averages over finite volumes (see also Section 5). The 
validity of the McGaugh (1991) calibration itself is uncertain. One reason for this lies in 
the ambiguity in the definition of U, discussed in the previous point. A second reason 
is that the calibration was made from a given set of ab initio simple photoionization 
models, and it is difficult to estimate a priori whether such calibration can be applied 
to real H ll regions. 



2.6.1. The Metallicity of NGC 5461 
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ci) A different metallicity calibrator could be used to resolve the degeneracy in R23, such 
as the [N 11] vs Ha indicator proposed by van Zee et al. (1998): 

12 + Log (0/H) = 1.02 Log ([N II] 6548, 6584/Ha) + 9.36, (6) 

which yields Log 0/H = 8.58 and 8.59 (that amount to Zg^s — 0.0105) with the Torres- 
Peimbert et al. (1989) and the Luridiana et al. (in preparation) data respectively; this 
result seem to confirm the hypothesis of a higher metallicity. 

We leave the metallicity question open for the moment, since a definite answer to it 
should come from different criteria than those taken into consideration in the present study. 
We will briefly discuss the variation of our models with metallicity in Section 6. 

2.7. Line Intensities 

The observed line intensities are listed in Table 3, together with other derived physical 
parameters. To compare with the models' predictions, we used the spectroscopic data by 
Garnett ct al. (1999), Torres- Peimbert et al. (1989), and Luridiana et al. (in preparation). 
Wc also used the /(H/?) value by Skillmann & Israel (1988) to map the nebular emission 
with larger apertures. 

Two of the slits considered (B and C) are very similar in size. Accordingly, we expect the 
corresponding data to be similar, at least as far as the slits have been placed on exactly the 
same position, and within the observational errors. Indeed, the line ratios /(A 5007)//(H/9) 
and /(A3727)//(H/3) are in good agreement; this is important, since these lines are the most 
intense, and they dominate the ionization structure of the nebula. /(A 6300)//(H/3) shows 
a difference of almost 50%, but this line represents a minor stage of ionization. If this line 
is formed in filaments and condensations, as it has been suggested (see, e.g., Stasihska & 
Schaerer 1999), one possible explanation for the difference between the data is that one slit 
captured more such condensations than the other. It should also be taken into account that 
the line is very weak. J(A 4074)//(H/9) and J(A4363)//(H/?), also weak lines, are discrepant, 
but they agree within the observational errors. 



11 



3. 



THE MODELING PROCEDURE 



3.1. 



The Numerical Models 



The photoionization models of NGC 5461 have been calculated with Cloudy 90 (ver- 
sion 90.05; Ferland 1996). We refer to the original documentation for a description of the 
characteristics of the code. The ionizing sources have been taken from the data set Star- 
burst99 (Leitherer et al. 1999), for the standard mass-loss case. 



Following the procedure outlined in Luridiana et al. (1999), we corrected the models' 
predictions for the slit size, before comparing them to the observational data. This is a 
necessary step, since a slit samples only a fraction of the nebula, in such a way that the 
ionization fractions, and thus the line intensity ratios, can be dramatically different from 
those of the complete model. The features of the slits used by Garnett et al. (1999) (slit 
A), Torres-Peimbert et al. (1989) (slit B), Luridiana et al. (in preparation) (slits C and D), 
and Skillmann & Israel (1988) (slit E), have been listed in Table 1. Given the small size of 
slit A, we expect the observational data by Garnett et al. (1999) to depend much more on 
the exact position and the small-scale structure of the nebula than the other data. Hence, 
in the comparison we will consider more significative a good fit with the other data sets. 



As already pointed out by Castaiieda et al. (1992), the modeling of a 3-D region from 
a 2-D image is an ill-defined problem, since many solutions are in principle possible. On 
the other hand, a totally realistic model of a region can never be calculated, because taking 
into account all the features of the small scale structure of a region rapidly outpowers any 
computational tool, not to mention the fact that observational data have a finite resolution 
and are affected by errors. 

To find a fair middle point between the simplistic and the nihilist points of view, one 
should ask himself which quantities are really relevant and worth (as well as possible) to 
determine, and which are not. With this criteria in mind, we adopted the following set of 
observational constraints: 

a) The relevant line intensity ratios. 



3.2. Comparison with the Observational Data 



4. 



THE OBSERVATIONAL CONSTRAINTS 
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b) EW(lip). 

c) /(H/5). 

d) The brightness profile (cfr. Castaiieda et al. 1992). 

e) The 7(A6717)//(A6731) ratio. 

f) The age range inferred by EW{}iP) (see Section 2.4). 

g) The degree of ionization. 

By 'relevant' line ratios we mean line ratios that fulfill as many as possible of the 
following requirements: i) The line is bright, so that the observational errors do not sensibly 
affect the analysis; ii) The atomic physics is well known; in) The line is produced by known 
mechanisms. By these criteria, the A 3727 and the A 5007 line intensities are very important 
constraints, while A 4363 is a less robust constraint since it almost surely has a contribution 
from processes other than photoionization. The A 6300 line is generally undcrprcdictcd in 
photoionization models, and the places and circumstances of its formation arc still being 
understood (see the discussion in Section 5.1.1). The [S ll] line intensities are important 
mappers of the low-ionization zone, but unfortunately the atomic physics of this ion is still 
not well known, and published parameters vary a lot. Finally, A 4686 is in principle an 
important constraint since it is a tracer of massive population, but it is subject to the many 
theoretical uncertainties still affecting our knowledge of W-R stars. EW{}iP) is an important 
constraint, because it is an age indicator, and it is measurable with a small error; however, 
its interpretation is limited by the circumstances mentioned in Section 2.4. /(H/3) is directly 
related to the ionizing power of the source, thus it is a constraint of primary importance. The 
brightness profile is a convolution of the ionizing radiation emitted and the gas distribution, 
so it is a basic constraint in any attempt to model the three-dimensional structure of a region. 
The [S ll] ratio is an important constraint, since it allows us to map the projected density of 
the region; the density dependence has to be deconvolved from the ionization-degree effects, 
giving clues to the low-energy range of the spectrum. 

Items a), b), c), d) and e) are quite strong constraints, since we tried to fulfill them as a 
function of the slit shape and size, so that each item actually consists of several interrelated 
constraints. Items c) and d) refer essentially to the same quantity, but averaged and displayed 
in different manners: in the first case, we aim at reproducing the observed fraction of /(H/3) 
intersected by each slit, while in the second case wc calculate the /(Ha) profile along a 
nebular diameter and compare it with the data by Castaneda et al. (1992). 
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5. RESULTS 

5.1. Best-fit Model Features 

Our best model has been calculated assuming a burst SF law of age 3.1 Myr, M^p — 80 
M0, and a Salpeter's IMF slope 1+x — 2.35. The rate of ionizing photons emitted, (5(H°), 
was set at the value (5(H°)= 3.47 x 10^^ photon s~^, as stated in Section 2.3. The stellar 
metaUicity has been set to — 0.008. 

The density law is given by A^e = 500 e'^^^"^^^^"^^ cm~^ in both halves, with r in parsecs, 
and 6 — ro — 54: pc (corresponding to 1.5" at the assumed distance of 7.4 Mpc). The filling 
factor has been set to e = 0.002 and e = 0.005 respectively, to reproduce the low and high 
brightness peaks seen in the Ha profile (Castaiieda et al. 1992). We stress again that a 
lower observed /(A6717)//(A6731) ratio does not necessarily imply a higher density, since 
the observed value depends on the contributions of all the gas parcels intercepted by a given 
line of sight, with weights depending on the local ionization and temperature conditions. 
The covering factor of the model is cf — 1. A central hole, of radius Tin = 20 pc and Tin = 35 
pc in the low- and high-e halves, is also present; the hole has only a minor effect on the line 
ratios, but it improves the fitting of the brightness profile between the peaks (see also the 
discussion in Section 2.2). The total radius takes different values in the two halves, averaging 
460 pc. The gas metallicity of the model has been set to Zg^s — 0.0066. Table 4 and Fig. 2 
show the features of this model. 

The first hne lists the logarithm of the total H/3 flux, in erg s~^, emitted by the complete 
model, and by the fractional volumes intercepted by the four slits. In the following lines we 
list the predicted line ratios, /(H/3), EW{}i(3), and R23, once again for the complete model 
and the four slit-biased cases. For each constraint, the boldface values in parentheses are 
the ratios between the computed and the observed values: i.e., a 'perfect' model should rate 
pure I's. 

5.1.1. Oxygen Lines 

The agreement is very good for A 5007 and A 3727, the most important lines according to 
the criteria set in Section 4. A 6300 is too weak in our model, a not unusual fact in the history 
of photoionization modeling, for which many possible explications have been invoked (e.g., 
Garcia- Vargas et al. 1997; Martin 1997; Stasihska & Leitherer 1996; Stasihska & Schaerer 
1999). 
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5.1.2. [S ii] A 6720 

A 6720 is predicted too weak by about a factor of about 2: more generally, in all 
our modeling attempts following an IBSF, we systematically found /(A 6720)//(A 3727) 
~ 0.5(/(A 6720)//(A 3727))°^^ whereas CSF models of ages ~ 5 Myr give /(A 6720) //(A 3727) 
- (/(A6720)//(A3727))°'^^ indicating that the IB spectra lack flux in the range 0.7 Ryd < 
i/ < 1.0 Ryd, range responsible for the production of S ll but not of O ll. The most prob- 
able explication for this result is the presence of an older and cooler stellar population in 
NGC 5461, not taken into account by the extreme IB scenario. 

5.1.3. i?ellA4686 

A 4686 is completely missing in our model, which has an age just prior to the W-R 
phase onset. The reason behind such a choice for the age is that the appearance of W-R 
stars in the IB scenario yields a sudden hardening of the spectrum, dramatically rising the 
ionization degree; furthermore, with any reasonable choice of the population parameters, 
even the A 4686 intensity rises too much. These circumstances might indicate that a strictly 
analytical treatment of the stellar population evolution is not realistic, and that the A 4686 
flux possibly comes from only one star, maybe slightly older than the average population. It 
should be also considered that the observed A 4686 flux also contains a stellar contribution: 
accounting for the stellar contribution, the observed nebular value is reduced by a factor of 
two. In the case of NGC 2363 (Luridiana et al. 1999), this effect was not taken into account 
because the stellar contribution is not important for the age and metallicity of that region. 
Summarizing, given the numerous uncertainties still affecting the theoretical modeling of W- 
R and W-R— like stars, the weakness of the observed A 4686 (see Table 3) and the additional 
uncertainties deriving from the statistical fluctuations expected in the high-mass tail of the 
mass distribution (Cervifio et al. 2000) A 4686 should not be considered a robust constraint 
for this object. 

5.1.1 [O III] A 4363 

The intensity of A 4363 with respect to H/? is slightly underpredicted by our model, with 
the exception of slit C (which carries the greatest observational error). The observational 
error on this line is 0.04 dex in the Torres- Peimbert et al. (1989) value, and twice that on 
the Luridiana et al. (in preparation) one, so that the two measurements agree within la. 
The average between the two observational values, taking the errors into account, would give 
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/(A4363)//(A4861) ~ 0.135, with a model/observation ratio of 0.12/0.135 ~ 0.88. However, 
if we consider the A 4363/ A 5007 ratio instead, we find that it is fairly well reproduced by 
the model. 

5.1.5. [S II] A 4074 

Finally, the weak [S ll] A 4074 hne is underpredicted by our model, but better reproduced 
than A 6720. 



5.1.6. Flux Fractions 

The line marked with /(H/3)''^'Y-^(H/3)*°* is calculated from the first line, and shows the 
fraction of the H/3 fiux intercepted by each sht, as compared to the total H/3 fiux emitted by 
the complete model. The values in parenthesis are the corresponding observational values, 
calculated as g(HO)/g(H°)'"^ = g(H°)/g(HO)S"tD. ^he agreement is very satisfactory. 



5.1.7. EW{}II3) 

The next line lists the predicted EW{Y{.I3) values. Again, the agreement is rather sat- 
isfactory, especially taking into account the uncertainties accompanying this quantity (see 
Section 2.4). 



5.1.8. R23 

Finally, the last line lists the calculated R23 values, which fit the observed values very 
well. i?23 turns out to be roughly constant for the different apertures, in agreement with the 
result found by Kennicutt & Garnett (1996) on the observational side. However, we must 
caution that this is true only as far as the observed volumes simultaneously sample low- 
and high-ionization zones, since the local R23 value spans along the nebula a range of more 
than one order of magnitude, reflecting the large variations of the local ionization parameter. 
Further resolved spectroscopic studies are needed to assess the line-of-sight variations in a 
real nebula. For the moment, it is safe to say that R23 should not be considered strictly 
constant, or, equivalently, that f/ is a local parameter. A commonly used definition of 
U, based on the inner conditions of the nebula, (i.e., U = Q{H^)/4:7:cRf^Nf,), is not very 
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useful, since it does not take into account the density distribution. As an extreme (but not 
unrealistic) example, when wc tunc R^^ from, say, 0.1 pc to 1 pc, the model stays essentially 
the same, while U varies by 2 orders of magnitude. 



5.1.9. Brightness and [S ll] emission profiles 

In Fig. 3 the predicted Ha and /(A6717)//(A6731) ratio profiles (solid lines) are com- 
pared to their observational counterparts. The zero-point of the computed profiles has been 
shifted to r = —1" (36 pc) to make it coincide with the zero-point set for the observational 
data in the original paper. 

The brightness profile was calculated assuming a sht width of 36 pc (1"). This is 
presumably the aperture used by Castaneda et al. (1992), as obtained by comparing their 
instrumental FWHM intensity (1.88 A) with the sulfur lines FWHM intensity (~ 2.8 A), 
and taking into account the dispersion (0.71 A pixel"^) and the spatial scale (0.33 " pixel^^) 
reported in the original paper. Fortunately, the profile is quite insensitive to the exact slit 
value used, at least for aperture values smaller than 3". We intentionally did not smooth 
the calculated profile, to show that the model faithfully reproduces the observed data in a 
very local sense; nevertheless, averaging the intensities values over 1" intervals, following the 
observational procedure, would undoubtedly improve the fitting. 

The agreement of our model with the observed /(A6717)//(A6731) values is also very 
good. The computed profile clearly shows how two exactly identical local density distribu- 
tions can give rise to different /(A6717)//(A6731) ratios, only by virtue of different e values. 
We also note that the observed /(A 6717)//(A 6731) value d — -1-4" is clearly unphysical. 



5.2. Total stelleir and ionized gas mass estimates 

The model described in Section 5.1 yields a stellar mass of Mf = 1.0 X 106 M© 

between 0.8 and 80 M©, or, equivalently, Mi^^^°^® = 0.9 x 10^ M© between 1.0 and 80 
M0. Rosa & Benvenuti (1994), observing NGC 5461 with a 1" aperture, fitted the observed 
continuum to that of population synthesis models, after correcting for the attenuation by dust 
and the nebular continuum emission. Following this method, they estimated a stellar mass of 
Mf >^-°^®(l") = 1.0 X 10^ M0 in the (2.0-80) Mq range, corresponding to m5>^-°^®(1") = 
1.4 X lO'"^ Mq with a Salpeter IMF slope. Adopting a scale factor of 4 to account for the small 
aperture used (Giannakopoulou-Chrcighton, Fich, & Wilson 1999), this figure translates into 
^m>i.oMq ^ g g ^ roughly 2/3 of our estimate. 
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Carigi, Colm, & Peimbert (1999) propose IMFs, based on the one by Kroupa, Trout, 
& Gilmore (1993), accounting for dark matter in the form of substellar bodies. Their IMFs 
are parameterized as a function of an r-value^, which depends on the assumed slope in the 
M < 0.5 M0 range. Their preferred IMF (corresponding to the r = 1.8 case: see their 
paper), truncated at Mup = 80 Mq, yields: 



0.206, 



(7) 



while with the IMF recently determined by Kroupa (2000) one obtains for the same mass 
range: 



0.377. 



These relations allow to estimate the total stellar mass of our model region: 



tot 



M>1.OM0 



/0.292 = 3.08 X lO^M 







(9) 



The total ionized mass obtained through straight integration of the local density over 
the volume is 

M*^ = 1.63 X lO^M©. (10) 

Observationally, this expression corresponds to the mass estimated through forbidden-hne 
density: 



M. 



FL 



A I eN^{FL)dv 
Jv 

Jv 



(11) 
(12) 

where A = rriH 4A^(He)/iV(H) + 1 Ar(H)/(^iV(H) + iV(Hc)j . A different observational 
estimate of the total mass of the ionized gas can be made by integrating the rms electronic 
density over the total volume: 



= ^ / Ne{rms)dV 
Jv 



(13) 



Peimbert (1966) showed that the expressions 12 and 13 represent lower and higher limits 
to the total ionized mass value, since in real nebulae the density contrast is not as extreme 



^This r has no relation with the r defined in Section 2.1 
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as supposed by the filling-factor scheme. In our case, this implies for the total mass of the 
ionized gas in NGC 5461: 

1.63 X lO^M© < Mg^° 5^61 < 3.25 x IO^Mq. (14) 

According to recent estimations by Giannakopoulou-Chreighton et al. (1999), the total 
molecular mass in NGC 5461 lies in the (15 — 40) x 10^ range, accompanied by 1-2 
times as much neutral mass. Taking an average value of 6 x 10^ Mq for the sum of these 
two components, and adding the ionized gas value, we find a total gaseous mass in the 
(6 — 9) X 10^ Mq range. The ratio of the total stellar mass to the total gaseous mass yields 
a star-formation efficiency in the (0.03-0.05) range, in agreement with current estimates of 
this parameter (e.g., Lada 1992; Evans & Lada 1991). 

We conclude by remarking that the rather high value found for M*°* imply that statistical 
fluctuations of the IMF do not play a significant role in this region (see also Cervino et al. 
2000). 

6. DISCUSSION 

In this Section, we wish to justify our choice for the parameters of our favoured model, 
by means of schematically illustrating the changes in the results obtained through variations 
in the input ingredients: 

6.1. Age 

The constraints on the age are set by the observed £'VF(H/5) value, and by the ionization 
degree of the nebula. For an IB scenario, the /(A5007)//(A3727) ratio steadily decreases 
with age until i ~ 3 Myr, with the exact age value depending on metallicity and on the 

other stellar population parameters. Then, the W-R stars are born and abruptly increase 
the ionization degree, which falls again at about 5 Myr following the death of W-R stars. 
Thus, in the — 5 Myr window there is only a short period around 3 Myr compatible with 
the observations. Ages greater than 5 Myr are excluded by the high EW{li(3) observed. 

The CSF case, which allows older ages, will be discussed in Section 6.5. 
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6.2. Metallicity 

The results presented here are scarcely dependent on metallicity, since at Z ~ Z0/3 the 
increase in the number of emitters is almost perfectly counterbalanced by the decrease in 
electron temperature, leading to almost constant oxygen line ratios. Different metallicities 
require slightly different age values, mainly because the evolution of the ionizing flux depends 
on mass-loss rates, which in turn depends on Z. 

6.3. Geometry 

The chosen geometry was the result of many crossed observational constraints: mainly 
the /(A 6717)//(A 6731) profile, e, the ionization degree, the nebula radius, and the brightness 
profile. Although we cannot ensure that the solution is unique, we are confident that at 
least qualitatively the real gas distribution is not too far from our model's assumptions. 
Nevertheless, it is useful to discuss possible variations in the input parameters determining 
geometry. 

6.3.1. Density 

A variation in the density normalization (iVo), leaving the remaining parameters un- 
changed, yields modifications in both the total radius and the ionization degree of the nebula. 
The total radius changes due to the constraint on (5(H°), the ionization degree changes due 
to the variations in the ionization and recombination rates. The recombination rate increases 
roughly with the squared density, while the ionization rate increases only linearly with the 
density, implying that a rise in A^"o at fixed radius yields a fall in the ionization degree. On 
the other hand, the ionization rate depends on the rate of photons striking a given point, 
which in turns depends also on the distance from the source; the average distance is smaller 
in a higher-density model if (5(H°) is fixed. In our configuration, the density effect outweighs 
the distance effect, so that an increase in A^o lowers the ionization degree. 

An interesting question that can be raised on this subject is how the model would change, 
had a simpler density structure been adopted (e.g., a constant density or hollow sphere, such 
as those available from grids of photoionization models) . The answer is that the features of 
our best model are strongly dependent on the density structure chosen. Given the size of the 
region, a constant density model must be characterized by a very low density, with the two 
opposite effects mentioned in the previous discussion (fall in the recombination rate due to 
the lower average density, fall in the ionization rate due to the higher average distance). As 
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a rule, for a given ionization source, the average ionization degree drops in constant density 
models. A second major change, with respect to our gaussian density-distribution models, 
is that the brightness profile is no longer reproduced. Several constant-density models were 
calculated to confirm this predictions. 

The implication is that, dropping both constraints on the region's radial structure (the 
/(A6717)//(A6731) and brightness profiles), a harder ionization source would be invoked to 
reproduce the observed fine intensities, leading to a strong bias in the inference of the central 
star cluster properties. Summarizing, we claim that to assess the properties of the ionizing 
field it is necessary to use tailored density distributions; this is perhaps the more important 
result of the present study. 

6.3.2. Filling Factor 

The filling factor acts on the ionization degree via the average distance from the ionizing 
stars. Since the recombination rate depends on the local (clump) density only, an increase in 
e, with the other parameters kept constant, lowering the average distance of the gas parcels 
from the ionizing source, yields a higher ionization degree. (Note that this is true only as 
far as the nebula is radiation bounded: the inverse trend can be found in density bounded 
objects with fixed radius, in which the higher the concentration towards the center, the less 
the low- ionization zone is 'sacrificed' by the constraint on the radius.) 

6.4. Density-bounded Models 

The density-bounded case can be described through the two limiting cases of a covering 
factor smaller than 1, and of a spherically symmetric nebula with R < Rs, where Rs is the 
Stromgren radius. 

In the first case the gas covers a solid angle Q < An (the covering factor is defined 
as cf — fl/in), i.e. the nebula is not spherically symmetric; the ionizing photons emitted 
in some directions are completely absorbed, while those emitted in other directions escape, 
resulting in a photon leakage independent of frequency. 

In the second case, the nebula is spherically symmetric, and the photon leakage af- 
fects preferentially the highest frequencies, i.e. those frequency characterizing the photons 
reaching further into the gas. The common features between these two situations are the 
following: a) the EW{}iP) is lower with respect to a radiation-bounded case, and b) if -F(H/3) 
is fixed as a constraint, (5(H°) must be increased in density-bounded models with respect to 
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ionization-bounded models. In the following, the two cases will be discussed separately; of 
course, real nebulae are intermediate between them. 

6.4-1- Covering Factor 

In the models with c/ < 1, we scaled the rate of ionizing photons according to the 
relationship (5'(H°) = (5°(H°)/c/, in order to preserve the total emitted H/3 flux. The 
consequences of assuming c/ < 1 in the models depend on which other parameters, if any, 
are correspondingly modified. This, in turns, depends on the obsevational constraints set. 
One of the most stringent constraints of the present work is the observed /(H/3) as a function 
of the aperture used. Starting from our reference model with cf — 1, any decrease in cf 
should be associated to a corresponding increase in Q(H°) to preserve the HjS flux seen 
through each slit. This would lead to an increase in the ionization degree, calling for some 
further change in the stellar source and/or in the gas geometry to be balanced. We found that 
moderate changes in cf (say, lowering cf from 1 to 0.75) do not alter sensibly the important 
line ratios. The only changes are the following: EW{}iP) decreases, getting farther from the 
observed values; /(H/3) slightly decreases, but the effect can be easily counterbalanced from 
a small modification in A^"™"; and /(A4074)//(A4861) improves due to the slight increase 
in temperature. On the other hand, /(A 6720) does not improve since the the ionization 
structure is essentially the same, and this line is less sensitive then A 4074 to changes in 
temperature. 

6-4-2- Eocternal Radius 

If the external radius is truncated before Rg is reached, the overall ionization degree of 
the nebula is altered, since the low-ionization lines are formed in the outskirts of the nebula. 
A detailed discussion on these models is quite complex, due to the interrelation of all the 
constraints, which implies that a change in one input parameter calls for changes in other 
parameters as well. Necessarily, our discussion will be simplified in this respect, and we 
will try to illustrate separately the cascade of consequences which results from relaxing the 
assumption R — Rg. 

The simplest possible solution is to simply truncate the nebula before the Stromgren 
radius, without further changes. The resulting models will generally preserve the value 
of most considered fine ratios, namely [O ill], [O ll], [S ll], and He ll fines. This is a 
consequence of the ionization structure of the nebula, of the density distribution, and of 
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the size of the slits as compared to the size of the region: these three factors contribute 
in the same direction (we are not considering here extreme cases of very small radii). The 
only line which is strongly affected is A 6300, which was already underpredicted. A major 
problem with models of this kind is that they do not reproduce the constraints on the 
observed brightness profile, which becomes too weak, and on the radius, which becomes too 
small. The constraint on the radius may then be fulfilled by a decrease in density and/or 
filling factor, but the problem with the observed brightness remains (it actually worsens); 
so, the only solution for truncated models appears to be to simultaneously increase Q(H°). 
But, again, the brightness profile is not reproduced if the matter distribution (A^e and e) 
is left unchanged; furthermore, e cannot be increased, since it would increase too much the 
degree of ionization. If is increased (by means of increasing either Nq or N^^'^), fairly 
satisfactory solutions can be found, as far as we do not depart too much from the reference 
radiation-bounded models. 

Summarizing, it appears that no substantial leakage of photons is affecting NGC 5461; 
as a gross estimate, we can say that at most 20% of photons escape from the region, and 
that the radiation-bounded model is a good approximation to the NGC 5461 case. Other 
solutions may still be possible, but they involve radical changes in the nature of the ionizing 
source (namely, a much softer radiation field), and we did not explore them in this work. 
And, of course, these results have been obtained for the particular case studied, while for 
other H ll regions the situation might well be different. 



6.5. SF Law 

As mentioned earlier, the observed EW{B.f3) values are compatible with both a young 
(t ~ 3 Myr) burst and a slightly older (t ~ 5 Myr) continuous SF event. A substantial 
improvement in the [S n]/[0 ll] ratio can be obtained with the CSF case, thanks to the 
contribution of older stars to the formation of [S ll] . However, in the CSF scenario the [O 
lll]/[0 ll] ratio remains too high for all ages and metallicity values, due to the continuous 
replenishment of hot, massive stars. We consider then that the IB is a better approximation 
to the real SF process going on in NGC 5461; an even better approximation would be obtained 
by adding to the IB spectrum the low-frequency emission from an older population. 
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6.6. IMF 

Since the slope of the IMF has no major effect on the results, we chose the standard 
Salpeter value 1 + x = 2.35. Regarding Mup, we chose a relatively low Mup value, since 
higher values yield far too high ionization degrees. For instance, if in our reference model we 
modify M^p from 80 Mq to 100 or to 120 M©, the [O iii]/[0 ll] ratio of the complete 
model changes from 0.42 to 1.37 or to 2.23 respectively; the trend followed by the slit-biased 
values is even more extreme, since, e.g., the [O in]/[0 ll] ratio seen through slit B changes 
in the same sequence from 1.14, to 5.64, to 13.41. We emphasize here that our results depend 
directly on the constraint enforced on the radial density distribution, and that many more 
solutions would be available if it were dropped. 

7. CONCLUSIONS 

We presented a selected photoionization model for the GEHR NGC 5461. The model 
is an asymmetric nebula, characterized by an off-center gaussian density distribution, with 
a peak value of = 500 cm^^ and different e values in the two components. The ionizing 
source is a young (3.1 Myr) burst with a Salpeter's IMF and M^p = 80 Mq, containing 4000 
O stars approximately corresponding to 3000 'equivalent 07 V stars, with the definition by 
Vacca (1994). The reasons underlying our choice of the IBSF reside mainly in the relatively 
low ionization degree of the nebula; CSF models are continuously replenished of hot, massive 
stars and maintain a high ionization degree. The age has been set taking into account the 
same constraint, and it is fully consistent with the observed EW{}iP) value. We estimate a 
total stellar mass of about M*°* = 3 x 10^ M© in the (0.01 — 80) M© range, and an ionized-gas 
mass lower limit of M*"* = 1.6 x 10^ Mq. Accounting for the gas in neutral and molecular 
form, we find that the star-formation efficiency lies in the 3 — 5% range. 

Our results are pretty robust with respect to reasonable variations in the input ingredi- 
ents. In particular, we are confident in our estimates of the stellar population's parameter, 
leaving open for the moment the questions of metallicity and W-R population. We are also 
confident that the overall matter distribution chosen (roughly, the amount of ionized gas 
found at each radius, resulting from the interplay of density and filling factor) is a good 
approximation to the real one, while we could not assess whether the constant-filling-factor 
configuration can be replaced by other gas distributions with different combinations of den- 
sity and fiUing factor. We consider that the region can be satisfactorily described as radiation 
bounded, and that no substantial leakage of photons is taking place. Stated synthetically, 
this comes as a consequence of the relatively low degree of ionization, coupled with the 
density and brightness profiles. 
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We successfully reproduced the following constraints: [O ill] /(A 5007) //(H/3) , A 4363/A 5007, 
and [O ll] /(A3727)//(H/3) as a function of the slit aperture; the [S ll] /(A 6717) //(A 6731) 
profile along the nebula; the Ha profile along the nebula; the equivalent width of H/5 as a 
function of the slit aperture; the R23 parameter. We failed to reproduce the observed He II 
A 4686, [O III] A 4363, [O l] A 6300, and [S ll] A 6720 line intensities. We attribute the failure 
in reproducing the [S ll] A 6720 line to an inadequate representation of the older, cooler part 
of the stellar population. This explanation could apply also to the A 4686 case: however, 
the low W-R statistics, and the theoretical problems still existing in the modehng of the 
atmosphere of these stars, might provide alternative explanations. 

Our study also demonstrates the following points: 

1. Photoionization modeling can be used to constrain the properties of H ll regions, as 
well as to assess our understanding of physical processes in ionized plasma. 

2. A detailed modeling should take into account as many constraints as possible, and the 
observed properties should be reproduced not only globally, but as locally as possible. 

3. Resolved spectroscopic studies are needed to provide modelers with large, homoge- 
neous, and exhaustive data sets for each modeled region. 

The authors acknowledge Grazyna Stasinska for many challenging comments, which 
greatly contributed to improve this work. V.L. would like to thank Daniel Schaerer and 
Corinne Charbonnel for extensive support, and the Observatoire Midi- Pyrenees de Toulouse 
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Table 2. Chemical abundances for the t'^ — case. 



Element 


Abundance^ 


N(He)/N(H) 


0.0904 


[0/H] 


-3.61 


[C/H] 


-3.78 


[N/H] 


-4.74 


[Ne/H] 


-4.29 


[S/H] 


-5.30 


[Ar/H] 


-5.82 




0.0066 



-[X,/R] = Log N{X,)/N{R)- 
Zgas is given by weight. 
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Table 3. Observational data. 



Observer 



Quantity 


Gal99^ 


TPPF89'' 


LEPOl^ 


LEPOl^ 


SI88'^ 




(Slit A) 


(Slit B) 


(Slit C) 


(Slit D) 


(Slit E) 


[O II] /(A3727)//(H/5) 


2.980 


2.138 


1.977 






[S II] /(A4074)//(H/?) 




0.013 


0.018 






[\J lllj I [A ^ouo ) / 1 [rip ) 


<U.U4U 


U.Uio 


U.Uii 






He a 7(A4686)//(H/3) 




0.021 








[O III] 7(A5007)/7(H/3) 


3.110 


3.020 


3.112 






[O i] /(A6300)//(H/3) 




0.015 


0.022 






[S II] /(A6720)/J(H/5) 


0.475 


0.225 




0.293 




[O III] /(A4363)/J(A5007) 


<0.013 


0.005 


0.004 






[S II] /(A 6717) //(A 6731) 




1.202 




1.321 




Log Q(HO) 




52.27 


52.18 




52.54 


EW(if/3) 




195 


157 






Log i?23° 


0.84 


0.79 


0.78 






Te([0 IIl])f 




9200 


8400 






Ne([S 11])^ 




370 




150 




Ne([0 11])^ 






190 







^Garnett et al. 1999. 

^Torres-Peimbert et al. 1989. 
'^Luridiana et al., in preparation. 

'^Skillmann & Israel 1988. 

ei?23= (/(AA 4959, 5007) + 7(A3727))/7(H/3). 

^Calculated by us with the observed line ratios and homogeneous atomic param- 
eters. 
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Table 4. Best-fit Model Predictions. 





Constraint 


Complete Model 


Slit A*^ 


Slit 


Slit C^ 


Slit D'^ 


Slit 


Log 


/(A4861) 


40.233 


38.496 


39.916 


39.935 


39.968 


40.162 


[O II 


/(A3727)//(H/3) 


3.679 


1.645 
(0.55) 


2.257 
(1.06) 


2.312 
(1.17) 






[S II 


7(A4074)/7(H/3) 


0.032 




0.011 
(0.91) 


0.011 
(0.64) 






[0 III] 


J(A4363)/7(H/3) 


0.007 


0.015 
(>0.37) 


0.012 
(0.79) 


0.012 
(1.05) 






[He II 


7(A4686)//(H/3) 


0.000 




0.000 
(0.00) 








[O III 


/(A5007)/7(H/3) 


1.553 


3.195 
(1.03) 


2.575 
(0.85) 


2.529 
(0.81) 






[Oi] 


J(A6300)//(H/3) 


0.023 




0.003 
(0.21) 


0.003 
(0.15) 






[S II 


/(A6720)//(H/?) 


0.449 


0.093 
(0.20) 


0.145 
(0.65) 




0.209 
(0.71) 




[O III 


/(A4363)//(A5007) 


0.005 


0.005 
(>0.38) 


0.005 
(1.00) 


0.005 
(1.20) 






[S II 


/(A6717)//(A6731) 


1.350 




1.228 
(1.02) 




1.280 
(0.97) 






^slit 1 ptot 






0.48 
(1.13) 


0.50 
(1.14) 




0.85 
(0.85 ) 




EW{11I3) 


263 
(1.06) 




126 
(0.65) 


132 
(0.84) 








-R23 




0.77 
(0.92 ) 


0.75 
(0.95) 


0.75 
(0.96) 







*Same aperture as used by Garnett et al. 1999. 
^Same aperture as used by Torres-Peimbert et al. 1989. 
'^Same aperture as used by Luridiana et al., in preparation. 
•^Same aperture as used by Skillmann & Israel 1988. 
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Fig. 1. — Density distribution of a typical model. A^o = 500 cm~^, 6 — 54: pc, vq — 54 pc. e 
is 0.002 on the left side, 0.005 on the right. The size and assumed position of the slits used 
in the subsequent analysis are also shown. 

Fig. 2. — Selected properties of our best-fit model. Upper panel: Ionization fractions of S+, 
S"*""^, and S^"*". Middle panel: [S ll] AA6717, 6731 emissivities as a function of radius; note 
how the density profile enhances the [S ll] contribution near the center, where the [S ll] 
ionization fraction is low. Lower panel: Density and filling factor as a function of radius. 

Fig. 3. — Upper panel: Ha brightness profile of the best-fit model (sohd line), com- 
pared to the observed profile by Castaneda et al. (1992) (dot-dashed fine). Lower panel: 
/(A6717)//(A6731) ratio profile of the best-fit model, superposed to the observational data 
by Castaneda et al. (1992) (asterisks). The zero-point of the computed profiles has been 
shifted to r = —1" to make it coincide with the zero-point of the observational data. 
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T 1 1 1 1 1 1 1 1 1 1 1 1 r 



Slit E 




Slit A: Garnett et al. (1999) 
Slit B: Torres-Peimbert et al.(1989) 
Slit C: Luridiana et al. (3000), blue 
Slit D: Luridiana et al. (2000), red 
Slit E: Skillmann fr Israel (1988) 
J I I I I I I I u I I I I I I 

-500 500 

Distance from center (parsecs) 
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Depth (parsecs) 
-200 200 




1 - 

I I I I I I I I I I I I I I I I I I I 
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Depth (arc sec) 



